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dimensions with different symmetries. In particular, we show that a U(l) LGT defined in a 
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is achieved by expressing the U(l) LGT partition function as an inner product between two 
quantum states. 
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1. Introduction 

The partition function 3f is the keystone in both statistical mechanics 03 and quantum field 
theory 0. Thermodynamic quantities, like the free energy, the entropy, correlation functions, 
etc. can be evaluated once is known as a function of physical parameters (such as the 
inverse temperature /3, couplings constants /, external fields and other order parameters). 
On the other hand, gauge theories have proven to be crucial in the description of nature. In 
particular, quantum electrodynamics (QED) is described by a U(l) gauge theory. 

It has recently been shown |f3l |U that the partition function of any classical spin model 
can be mapped to that of an (enlarged) four-dimensional (AD) lattice gauge theory with gauge 
group Z2 (see also the original idea and closely related works (6l|7l[8l). More precisely, if 
one tunes the coupling strengths of the partition function of a (large enough) AD Z2 LGT, this 
would equal the partition function of a classical spin model in any dimension, with any type 
of interaction pattern (including arbitrary many-body interactions), and thus also includes 
models with local and global symmetries. To obtain this, one expresses the partition function 
of a large class of models as an inner product between two quantum states, and then relates 
the quantum states (see JU for a detailed treatment of the quantum formulation). The result 
is very general and unifies very different models. In this sense, the AD Z2 LGT is a complete 
model for this class of discrete models. In fact, these mappings from a general discrete model 
to a Z2 LGT allow one to gain structural insight. 

In the present work, we show that the partition function of many continuous classical 
models can be expressed approximately (to arbitrary precision) as the partition function of the 
AD U(l) LGT. The result holds exactly for a large class of models, namely, models whose 
Hamiltonian has a finite Fourier series and no constraints on the variables. This includes 
models of different geometry, and in arbitrary dimensions (ID, 2D, 3D, etc.) In the proof 
of the statement we generalize the quantum formulation developed previously j6l HI [TUl El 
and generate the (truncated) Fourier series of any target Hamiltonian. In parallel with the 
AD Z2 LGT as the complete model for discrete models, the AD U(l) LGT in this sense is 
complete for a certain class of continuous models. We also define and consider the U(l) 
LGT defined in a curved spacetime background. We show that this class of models are also 
included in our completeness result, i.e. they can also be mapped to the flat AD U(l) LGT. 
The physical content of this result is the following: as long as the gravity coupled to the QED 
is a fixed background, it can be absorbed as if it were a flat spacetime as far as completeness 
is concerned. This situation has a physical counterpart in cosmology where one finds photons 
in an approximate fixed curved spacetime. 

In this paper we will first present some basics on the U(l) LGT in §|2L Then we will 
prove the completeness of the U(l) LGT in §j3j. Further illustrations of the completeness 
result with some examples and applications will be given in §0]. We will proceed to the U ( 1 ) 
LGT in a curved spacetime and relate to the main result in §[5]. In §[6] we will generalize our 
result to a larger class of models. Finally, the conclusions will be drawn in §[7]. 
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2. Basics on the U(l) Lattice Gauge Theory 

LGTs are gauge theories on a lattice representing a discrete spacetime. Generically, LGTs 
are useful non-perturbative formulations of gauge theories which allow for numerical 
simulations, e.g. using the Monte Carlo methods |fTTfl. This allows one to go beyond 
perturbative calculations with Feynman diagrams. The abelian U{\) LGT was introduced 
by Wilson |[T2l| and Polyakov lfT3l [T4| as a generalization of Wegner's Ising gauge theories 
lfT5l . The U{\) LGT can be considered a discretization of electrodynamics (a pure gauge 
theory with an abelian gauge group U(l)) defined on discrete spacetime lfT6l . Non-abelian 
continuous LGTs are successful as having been proven to be asymptotically free [fT71 in the 
weak coupling limit, and important in the study of quark confinement at the strong coupling 
limit lfT2ll . 

For abelian LGTs it is possible to simulate numerically the continuum limit, A — > 0, in an 
appropriate way, and verify that the resulting theory is the actual QED without confinement. 
This verification is important since on the lattice, most LGTs show confinement which is an 
artifact. It so happens that in the strong coupling limit of an LGT, the property that the gauge 
group is compact is essential for observing confinement, regardless of whether it is abelian or 
non-abelian. Related to this, the phase diagram of LGTs is very rich and relevant for taking 
the continuum limit [16], and the U(l) LGT serves as a test ground for this purpose. Unlike 
the non-abelian SU(N) LGTs, which have a discrete center subgroup Z#, the abelian U(l) 
LGT has a continuous center subgroup which is identical to the group itself. The role of the 
t/(l) group on the confinement to deconfinement transition is also of interest from the point of 
view of the abelian dominance hypothesis for confinement which holds that a U(l) subgroup 
controls the non-perturbative dynamics of non-abelian gauge theories [fT8l . 

In the mean time LGTs have emerged as interesting theories by themselves. They are 
examples of models with local symmetries and non-local order parameters. They exhibit 
phases which do not appear in the continuum limit lfT6ll . The phases of the U(l) LGT can be 
characterized by the Wilson loop (see £0 below) which is an gauge invariant order parameter 
of the model. In the confined phase, this order parameter obeys an area law, whereas in the 
unconfined phase it obeys a perimeter law. The 't Hooft loop is a dual variable to the Wilson 
loop, and it constitutes another order parameter. 

There are compact and non-compact 4D U(l) LGTs QTJ [TH |20l |2T||. In the compact 
U(l) LGT the degrees of freedom are exponentials of the edge degrees of freedom. Thus, 
they are directly elements of the U(l) group which is compact, thereby the name. On the 
other hand in the non-compact U(l) LGT the degrees of freedom are associated directly to 
the edges of the lattice, which do not need to be elements of the U(l) group. The former 
one only has a deconfined phase with massless photons, thus one recovers QED in the limit 
of continuous spacetime. The latter one has two phases: a weak coupling phase, where the 
model has massless photons (gapless excitations), and a strong coupling phase, where there 
are massive photons and magnetic monopoles. The photons are screened by the monopoles, 
which corresponds to a mechanism of confinement of electrical charge. These phases are 
characterized with Wilson loops and 't Hooft loops as mentioned above. In this paper, we will 
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deal with the compact U ( 1 ) LGT. 

We briefly summarize the formulation of the U(l) LGT as follows (with our notation 
close to that of lfT6lO . It is illustrative to derive the Lagrangian of classical electrodynamics 
from that of U(l) LGT in the limit of continuous spacetime ifToll . Let us consider a AD square 
lattice as the discretized spacetime. Let A denote the lattice spacing and a, b etc. the unit basis 
vectors. We denote vertices, edges and faces of the lattice by v, e and /, respectively. The set 
of all vertices, edges and faces is denoted by V , E, and F and the number of elements in each 
set by |V|, \E\, and |F| accordingly. We choose a direction for each edge, see Fig. [Q The 



P + Ab P+Aa+Ab 



n+b n+a+b 




hi v 



P a P + Aa 
(a) 




Figure 1. (a) A face of the square lattice, a and b are unit basis vectors and each edge is 
assigned a direction ±a or +b. The four vertices of the face are labeled by P, P + aA, P + bA 
and P + aA + bA. (b) For notational convenience, we use a single labeling for the vertices, 
namely, n (the « th vertex), n + a,n + b and n + a + b, with a = aA and b = bA. 

gauge field A a (n) is defined along each edge incident to vertex n with direction a indicated by 
the sub-index (see Fig. |2]). The Wilson loop for an elementary face (such as in Fig. HJa)) is 

U = c iA ^M^ c ^M^+^) e -^MP+^b) e -iM A h {P) 

= e iA^([A b {P+M)-A b (P)]-[A a {P+Ab)-A a (P)]) ^ 

in which ^ is the coupling constant. The exponential of the second line in (0Q) resembles 
the field tensor F ab = d a A b — dbA a . The action is constructed from the Wilson loops, for 
notational convenience, we switch to the simplified notation as in Fig. UXb). We make the 
following redefinition of the gauge field: 

6 a (n)=Am a (P), (2) 

as well as the following conventions: 

6- a (n + a) = - G a (n), 

& ab (n) =e a (n) + e b (n + a) + 0- a (n + a + b) + 0^ h (n + b). (3) 

We can take the redefined gauge field a as defined on the edges instead of on the vertices, 
see Fig. [2l With these redefined gauge field notations, the exponential in (OQ) is equal to i& ab 
(which is a sum over the four 6 e field variables along the edges of a face, e E f). & ab can be 
thought of as defined for each face and denoted alternatively by : 

0/ = 0i + 2 - 03 " % = £ XeO e . (4) 

eedf 
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6_ a {n + a + b) 



6_ b (n+b) V 



(a) 



A e b {n + a) 




0, 
(b) 



Figure 2. (a) The gauge fields Q a as denned along the edges, (b) A further simplified notation 
where the fields along the four edges of a face are denoted by 0\, O2, 63 and 64, and Q e in 
general. 



where df denotes the boundary of face /. The signs of the four edge variables in © 
are chosen to be consistent with Fig. [2] (b). In general, a face can contain edges pointing 
clockwise and others counterclockwise. So that % e = 1 [—1] if e is oriented clockwise[counter- 
clockwise]. (Note that only the relative sign of the e 's is relevant for the action © below.) 
Now the Wilson-Kogut action of the U(l) LGT can be written as 

5 =il [l-cose/]. (5) 

^ f€F 

This (Euclidean) action in the naive continuum limit is seen to be £ f eF ^ ~ \ J d 4 x F^ b , 
consistent with the classical 2/(1) theory. Also, the action takes the same form as in © in any 
spacetime dimension. This action is invariant under the gauge rotation applied on any vertex 
v: 

gv= n u ? e i ^ e ( 6 ) 

e: v(zde 

where U e is an element of the gauge group, in our case, U e G U(l), and de denotes the 
boundary of e (i.e. the product in © applies to all edges e incident to vertex v). Finally, the 
partition function (which corresponds to the Euclidean path integral) of this pure U(l) gauge 
theory on a lattice takes the following form: 



(7) 

K e<EE ) y " J feF 

Note that the first constant term in the sum of the exponential in © only introduces an 
overall constant factor which can be omitted^ In addition, the coupling constant ^ can be 
made face-dependent. Accordingly, we will start with the partition function of the following 
inhomogeneous model with local coupling constant Jf for each face (The coupling constant 

I Ignoring constant factors of the partition function corresponds to an overall shift in the action (or the 
Hamiltonian) which never changes any observable nor the equations of motion. It only introduces an additional 
constant to the free energy as we take the logarithm. Observables only relate to differences (or derivatives) in 
free energy. 



U(l) LGT Connects All Classical Models with Continuous Variables 



6 



Sf and (possibly) an inverse temperature /3 are all contained in the local //'s): 




3. 4D U(l) Lattice Gauge Theory as a Complete Model 

In this section we show that the AD U{\) LGT is complete in the sense that a large class of 
classical partition functions with continuous degrees of freedom can be represented as special 
cases of the partition function of this model. In § 13.11 we define the set of models considered 
in our completeness result. Then, we present some tools in § 13.21 and § 13.31 that we will require 
for the proof of the main result. Finally, our completeness result is proved in § 13.41 . 



3.1. Class of Classical Models 

Our completeness results embrace all continuous classical 'spin' models, i.e. models with a 
Hamiltonian satisfying the following conditions: 

(i) The Hamiltonian depends on a set of N continuous real variables {xj\j — 1,2,. ..,N} and 
each variable takes value within a finite interval, i.e. 

Xje[aj,bj], Vj. (9) 

(ii) The Hamiltonian is a sum over ^-body interactions with 1 < K < N, i.e. 

■*({*;})=£ I " w ({*./})■ do) 

K=l{K-body} 

(iii) There is no constraint on this set of variables, i.e. 

Xj with j = 1,2,..., N are independent variables. (1 1) 

Furthermore, we assume that each ^-body Hamiltonian is a well-behaved function which 
allows a Fourier series expansion over the Xj variables. (Note that the spacial dimension of 
models under consideration is arbitrary.) 

With the above conditions, we first normalize the ranges of variables by a linear change 
of variables from {xj} to {6j} 

1% ( a; + bj\ 
0j = 7—lxj--^-) i Vj (12) 



b j~ a j 

such that each 6j G [— %, %] . The relation (fT2l) can be inverted Xj = bj 2 ^ J 9j + aj t bj so that the 
Hamiltonian is reexpressed in terms of these normalized variables and denoted by J4? ({0/}) • 
Next, we make a Fourier series expansion of each if-body Hamiltonian ({0j}) over the 
6j variables: 

/ K terms 

({Gj})= £H« } exp i ZmjBj 

{>nj} \ j 

j K terms \ / K terms \ 

= I ^ H {fn!} C0S I m J & J -SH^sin £ mj0j .(13) 

{mj} \ j ) \ j / 
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where mj G Z, Vj and the Fourier coefficients 9tH| m ^ and — 3HV are all real (91 and 3 
denotes real and imaginary parts, respectively). Therefore the general Hamiltonian Jt? in (fTDl) 
is now written as a Fourier series over a set of basis functions 

\mj\ = 0,1,2,..., ;=1,. ..,#}. (14) 



{cos j 




| , sin | 













5.2. Quantum Formulation 

We present a quantum representation of the partition function 3f in ©. First, assign a 
quantum state |0y) to each face: 

I®/> = IL^)- (15) 
eef 

Then we define the following quantum state | \jf) which contains the interaction pattern (By 
'interaction pattern' we mean the lattice or graph structure representing which variables are 
interacting): 

iv>= f* (n de <)®\ & f)- < i6 > 

J ~* \eeE J f eF 

Next we define another state | a) which contains the (Euclidean) weight for each configuration 
of the e variables of the whole lattice: 




\a)= I [ [[ dd e I I ll e/C0 f )Q9\®f)- ( 17 ) 

/ feF 

In the following we show that (oc\\j/) is proportional to the partition function iF in ([8]). 
Let us consider the linear transformation from the edge variables (gauge fields) 9 e to the face 
variables (field tensors) 0^. We define two column vectors: 

0= (0 l ...0 f ...0\ F \) t . (18) 

where ' denotes transposition. They are related by a linear transformation = 1-0 where I is 
the face-edge incidence matrix with matrix elements If >e = 1 if e e df and otherwise. 

We always assume that all edge variables are independent, i.e. rankl = For a square 
lattice with periodic boundary conditions in D-dimension (D > 2), there is a relation between 
\E\ and |F| such that |F| = ^^-\E\. If\F\ < \E\ we add a set of linearly independent auxiliary 
face variables F aux , where |F aux | = l^l — This results in an l^l x \E\ incidence matrix I, 
viz. 

0= (0i... 0| F |...0| £ |) f = I-e. (19) 

Geometrically, this amounts to saying that we are introducing new auxiliary faces to the 
original lattice and these new faces need not be squares. If \F\ > \E\, we pick up a maximally 
independent set of face variables and eliminate linearly dependent rows until the incidence 
matrix has a size l^l x So that we may write = (©i . . -0|£|) f =1-0, with again an 
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\E\ x \E\ incidence matrix I. Because not all face variables appearing in the quantum states 
| iff) (fT6l) and | a) (fTTT) are independent, when we take the inner product of | iff) with | a) , formal 
infinities would arise. These infinities, e.g. 8 2 (&f — 0^-) = 5(0/ — 0^)5(0), always take the 
form of powers of 5(0). This problem can be treated in the following way. We introduce a 
large 'momentum' truncation A in the Fourier space (inverse space) for each face variable. So 
that 5(0 = 0) = 2^ LrceZ e1 " ~ 2^ with A — > °°. Therefore when taking the inner product we 

\F I \E\ 

shall have an overall factor (^) which tends to infinity formally. However, this factor 

does not affect any observable because it only introduces an additive constant to the free 
energy as we take the logarithm of the partition function. We shall call this a 'regularization' 
method that allows us to make sense of the expressions that come out with harmless infinities. 
(A similar constant factor has been omitted going from © to ([8]), see also (1251) below and 
the argument there for comparison.) Finally we see that (a\ yf) is proportional to the partition 
function J° in (ED): 



(a\\ff) 



Udo' e do e ) n e 

~ n \ eeE J \feF 



1 r% 



detl 



I-n (n. 

(A) 



J feF 



feF am d& f 



1 f A N 

I detl I 



3T, 



3? 



detl v 1 



E\-\F\ 



\F\ < \E\ 



\F\ > \E\ 



const, x 3f . 



(20) 



3.3. Construction of the Fourier Basis Functions 

Here we show how to generate Fourier series expansion (fT3l) from a AD U(l) LGT, that is, 
how to generate the basis functions of (fl4l) . In this way a subsystem of our complete model 
will behave as the target model; more precisely, the Hamiltonian of the complete model on 
this subsystem will coincide with the (Fourier series of the) Hamiltonian of the target model. 
We remark, though, that we will generate the basis functions of (fl4l) only for bounded values 
of rrij, < mj < Mj for all j, where the Mj's are large, positive integers. We will show this by 
first introducing the merge and deletion rules, and by then explaining how to obtain arbitrary 
many-body interactions and build the Fourier basis from these fundamental rules (an explicit 



construction of many-body interactions with more technical details is given in Appendix A). 

Merge and deletion rules. For simplicity, we consider the state \a) of (fT71) defined only 
on two faces (see Fig.[3ta)): 



\u a ,b) = J_ x [U dd ej e y « cos0 «e^ cos0 "|0 fl )|0 fc ) . (21) 
We define the merge rule on, say, face a. Our aim is to obtain a delta function 5(0 a ), in the 
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limit J a — > oo. We consider a slight modification of (1211) : 

\6c a ,b)=N(J a )j 7t ( fl d0^® fl )e / -( cos8 »- 1 )ef'» CO88 *|® fl )|® 6 ), (22) 



where iV(/ a ) is defined as 



This implies that we will obtain the target partition function with the prefactors: 

terms „J m \ 

where M is the number of faces where the merge rule has been applied, and J m denotes the 
coupling strength of the merged face, J m — > °°. The usual quantity of interest, such as the 
free energy per particle in the thermodynamic limit, is shifted by a known amount (again, 
thermodynamic quantities involving derivatives of the free energy are not altered using our 
regularization mefhodj|): 

l im ZL = lim I i n z + Urn - ( Yj m - - In l^lA (25) 

where iV is the number of particles in the classical U(l) LGT model. 
We have 



N(J a ) f* d® a e J ^ & »-^ = v^I (J a ), (26) 

J-K Ja 

where /o(/ a ) is the modified Bessel function of the first kind. The asymptotic behavior of this 
function is 

/ (/ a )~ -^={l + ff{y)), Ja^°°- (27) 
V IK J a J a 

Thus it follows that 

lim N(J a )e Ja( > C0S@a -V = 8(® a ) , (28) 
as desired. Therefore, 

lim I" I FT de e )N(J a )e J ^ cos& ''-^e JbCOS&b \® a )\® b ) 

= I ddiddidd A de 5 de 6 de 1 & JbCO *® b \o)\® b ) . (29) 

J — 71 

where the variable ® b has become ® b after imposing the constraint of the delta function: 

0^ = 05 + 06 + 07-01-04-03- (30) 

§ The prefactor of the partition function in d24l . as well as the additive extra term to the free energy in ( l25l ). 
are all formally infinities. As far as these infinities appear in a controllable way, we can get them off by first 
choosing the coupling strength /„,'s large but still finite. So that a shift in the free energy has been introduced 
such as in (l25l l. Then we can set 7,,,'s to infinity after calculation of observables. No observable will be affected 
because they all involve a difference or derivative in the free energy. Since the additional term does not contain 
any physical parameter, it will not contribute even being formally infinite. 
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(a) 



05 







j — '—>» 

< 

<~m 


* < 
< 

•A34567 

-=-► 



<»0 6 



J = oo 





(b) 




(c) 

Figure 3. A blue square represents an interaction in that face, whereas a shaded blue face 
represents an face with an infinite coupling strength, (a) Merge rule: by letting Jf — > °° of the 
left face, this is merged with the face on its right. The resulting face only depends on the spins 
on its boundary, (b) The resulting face can be merged again with a neighboring face by letting 
Jf — >■ oo. Note the resulting direction of the arrows, (c) Deletion rule: (the interaction on) a 
face is deleted by setting Jf = 0. 



Note that this corresponds to a 6-body interaction of the same type, and with interaction 
strength J b . Thus, letting J a — > °° effectively merges face a and b. Note also that the variables 
with opposite sign have opposite pointing directions in the resulting face, see Fig.Oa). 

It is straightforward to see that the same derivation applies for \a) defined on all faces. 
In this case, one would only substitute the condition of the delta function, & a = 0, on the face 
with which it has to be merged. The process can be concatenated, that is, the face resulting 
from a merge rule can be merged again with a neighboring face, thereby becoming a 8-body 
interaction (see Fig. Ob)). Note that, whenever one face is merged to another, the new edges 
have the opposite direction than the original ones, as noted above. 

The deletion rule works by setting the Jf = 0, which results in switching off the 
interaction in that face (see Fig. [3£c)). 

Construction of many-body interactions. The Fourier basis functions of the set (fl4l) can 
be generated by making repeated use of the fundamental merge and deletion rules presented 
above. More specifically, one first generates several "copies" of a certain variable, that is, one 
replicates, say, mj times the variable 8j, as required in (fl4l) . Then one generates the cosine 
and the sine of an arbitrary sum of such term (also with the corresponding signs). We refer the 



reader to Appendix A for technical details of these explicit constructions. Here we emphasize 



that this construction is achieved by applying the merge and the deletion rule on specific faces 
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of the AD lattice, and by fixing some variables using the gauge symmetry of the model ("gauge 
fixing"). The latter is a procedure that can be carried out so long as the edges whose variable 
has been fixed do not form a closed loop ||2~2~1|. This is precisely why we need to resort to a 
AD lattice, since only then is our construction of the interactions (fT4l) free of closed loops (the 
same case as in j3l). 



3.4. The Completeness Result 

We are now in the position of proving the main result of this paper. In § 13.21 we have expressed 
in general the U{\) LGT partition function as an overlap between two quantum states. In § 13.31 
we have shown that through merge, deletion and gauge fixing in a AD U(l) LGT, any basis 



functions of a Fourier series expansion can be generated (see also Appendix A). As a result, 
any partition function of the U ( 1 ) LGT (with local coupling constants) of the following form 
can be generated and expressed in terms of a quantum amplitude (oc\\j/): 

^LGT (31) 

= f K (U 66 ^ eX P | L } 7 m; COS (X>;V) +4; sm [% m j Q J 
with 

0<mj<Mj, V/. (32) 

Here all Mj are (large) positive integers. Furthermore, in § 13.11 we have also expressed 
a general Hamiltonian with continuous variables in terms of a Fourier series with real 
coefficients. Let us compare (|3T| ) with the partition function of a general classical model 
satisfying conditions ©, (flOl) and (fTTI) : 

•^classical (33) 



ex p| E -/^ h k} cos (ll m j e ^j +/ 33H { m/ } sin ^L m / ^ 

with jV = Y[j 3 2% 3 • ^ e realize that if we choose the local coupling constants in (I3TT) to be 
equal to the Fourier coefficients (times /3), i.e. 

J c m , = -j89tH {m . } , J s m . = j83H {m . } , (34) 

then (up to a constant factor ,jV) the general partition function ^classical can be approximated 
by the partition function ^lgt- The approximation lies in the fact that in (1331) the integers nij 
run from to +°°, while in (|3Tb it runs within a finite range, see (|32l . This approximation 
can be made to any precision as we increase the repetition of the 6j variables in the basis 
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functions in (fl4l) . see | Appendix B[ Therefore, the AD U(l) LGT is complete such that the 



partition function of any classical model of continuous variables (without constraints) can be 
mapped to that of the former. That is, by fixing some of the coupling strengthes of a AD U{\) 
LGT, the remaining subsystem behaves as the continuous model to be simulated. 

We emphasize here that models satisfying conditions ©, (flOl) and (ITT]) are all included 
in the completeness result, regardless of their dimensions (can be larger, equal to, or smaller 
than D = 4) and specific forms of interaction. Also, though the original AD U(l) LGT is 
gauge invariant ©, this 2/(1) gauge symmetry is broken after the mapping because the target 
model no longer possesses the same gauge symmetry in general. 

Finally, we discuss the overhead in the system size of the complete model as a function 
of the features of the target model. That is, we study how many variables are needed in the 
complete model in order to generate the Fourier series of (TUT) with a truncation at the M th 



mode (each m 7 < M in (I321). We have seen in Appendix A that the generation of each Fourier 



component (each element of (1141) ) requires a polynomial enlargement of our complete model. 
i.e. In order to generate a single Fourier term such as cos(£ jinjOj), a number of Poly(Y,j m j) 
edge variables are needed, where Poly(-) is a polynomial of its argument. (This fact also 
implies that the same order Poly(Y,j m j) of couplings need to be tuned to infinity or zero in 
producing this Fourier term. See § 14.11 for a simple and explicit example.) Thus, the efficiency 
is measured by the number of Fourier components to be generated in the expansion d3jj). For 
a single ^-body interaction term H K ({0j}) in the Hamiltonian (depending on a certain set 
of K variables {$/}), we need to generate ~ [Poly(M)] K number of Fourier components, 
where Poly(M) is a polynomial in M. Generally, all combinations of Oj variables in K- 
body interactions, with < K < N (N is the total number of variables), may be present, thus 
resulting in the scaling 



[Poly(M)f ~ exp (N) [Poly(M)f n ™ (35) 



for the number of Fourier components. However, in most cases K does not scale with the 
system size, K = K max , and moreover only few-body interactions appear in the Hamiltonian 
(e.g. K = 2 for two body interactions). In this case the scaling is polynomial in both 
parameters 

? ) [Poly{M)] K ™* ~ Poly(N) [Poly(M)] Km ™ . (36) 



Note that this question of efficiency in generating Fourier components is unrelated to the 



question of the accuracy of the approximation of a finite Fourier series (see Appendix B for 
the latter). 



4. Examples and Applications of the Completeness Result 



As an illustration of our completeness result, we give a few example of models whose partition 
function can be reduced to our class of models. 
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4.1. TheXY model 

In the XY model we have 2D unit vectors (classical spins) ? ; defined on a lattice. The 
Hamiltonian obeys 0(2) (or U(l)) symmetry: 

J#xy = - £ Jij Si ■ sj = - £ 7 i7 cos (0 ; - Bj) , (37) 
<!';> <0> 
The partition function of the XY model 

2? XY = (n d£ ^ e^<u>JiM8i-0j) (38) 

is seen by itself to be of a U(l) LGT type. Therefore, for this particular case, we do not need 
to generate the whole Fourier basis, but we just need to generate the 2-body interactions of 
(1371) with the constructions explained in § 13.31 (more details in | Appendix A[). Here below 



we give an explicit pictorial construction of a ID XY model from a AD U(l) LGT as a 
transparent example of our completeness result. In this particularly simple case we only 
need a 3D sublattice of the AD U(l) LGT for the construction, as shown in Fig. H] below. 
The thick black edges represent variables of the target XY model which eventually will build 
a ID chain (in terms of interactions). All the odd numbered variables 9\, 83, %, etc. are 
distributed along the same line. Even numbered ones O2, O4, etc. axe, distributed separately 
along the two sides. This arrangement of variables guarantees the correct relative sign in the 
interaction cos(0, — (to be generated after merge of faces) between every neighboring 
pair of variables. In the figure it is illustrated how to obtain interactions cos($i — O2) and 
cos (02 — 63) by merging the blue faces and gauge fixing the red edge variables, and these are 
the only two typical constructions we need. By direct repetitions all interaction terms of the 
XY Hamiltonian can be constructed. It is easy to see that in producing an interaction term like 
cos(0j — 62), 12 couplings are taken to the infinity limit; similarly in producing cos($2 — O3), 
8 couplings are taken to infinity, etc. 

4.2. The Gaussian and Mean Spherical Models 

In both the Gaussian and the mean spherical models the Hamiltonian (also the 'effective' 
Hamiltonian appearing in the expression of partition functions) is a quadratic form of 
continuous unbounded spin variables (i = 1,2,... ,N): 

■^ = -% Jij<*i°j-Y,hi<*i- (39) 

<ij> i 

In the Gaussian model, each variable a,- takes values in (— +00) and is assigned a 
probability distribution of a Gaussian form: 

p{Oi)Aoi = J 2 q- Ao > da,-, i = 1,2, . . . ,N, (40) 
so that the average value is (of) — The partition function of the Gaussian model reads 

N 



(41) 
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Figure 4. Generation of the ID XY model from a 4D U(l) LGT. Red edges indicate edges 
whose variables have been gauge fixed, thick black edges indicate edges containing variables 
that are present in the target model (they will build the ID chain). Black arrow indicate cubes 
where the variable is replicated according to Fig. IA1 1 (Appendix A I. Interactions of the kind 
cos(0, — take place in blue prisms, and are a particular case of the interaction shown 



in Fig. IA4I ( Appendix A i. Indication of merged faces (blue faces in the figures) has been 



simplified here to avoid overloading. 



Alternatively, in the mean spherical model, instead of introducing a probability 
distribution, we impose a constraint on the average value of of such that 



N 



(42) 



(=l 



This constraint can be incorporated in the Hamiltonian by introducing a spherical field A : 

^iean = " £ JijOiOj ~ J>Oj + A £ of . (43) 

<ij> i i 

The partition function of the mean spherical model is now 



subject to the constraint 

<?^mean v 
\ 51 / 



1 dlnXv 



N. 



(44) 



(45) 



dX 1 /3 dX 

Comparing the form (l44l) with (|4TI) . we say that the Gaussian model is a mean spherical model 
with a prescribed spherical field A. 

For both the Gaussian and the mean spherical models, we can make a cut-off of the 
Oi variables in order to satisfy our condition ©. That is to say, we assign a joint probability 
distribution p{(7j} with a compact support ^ to the set of a,- variables (any Oj vanishes outside 
some (large) radius in M. N ): 

p{a,} = 0, VCii^, Mi. (46) 

As a result the integrals over <j/'s in (|4TI) and (|44b do not extend to infinity: 



(47) 
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With this modification of the models, both the Gaussian and the mean spherical models belong 
to our class which can be approximated by a AD U(l) LGT (the quadratic forms in (|4TT) and 
(|44|) all have Fourier series expansions). 

4.3. A Model of Coupled Planar Pendulums 

Consider a number of planar pendulums with pairwise coupling (a discretization of the sine- 
Gordon model in 1 + 1 dimensions), see Fig. [51 Let us assume that the interaction potential 




Figure 5. Coupled planar pendulums. The position (x,,;y,) of each pendulum is determined 
by the center coordinate (x,o,y,o) and an angle 0,. Pendulums are interacting pairwise with the 
interaction potential depending on the relative distances between their positions. 

between any pair of pendulums depends only on their relative distance and each pendulum 
has a self-interaction depending only on its position (e.g. a constant gravitational field). The 
Hamiltonian can be expressed as 

J$? = ^ w;(cos Of, sin 6j) + ^ u ij ( cos ®U cos ®ji sm ®U sm 0/> cos(0 ; - — Oj) X48) 

i <ij> 

where in the simplest case the self-interaction 

«/oe Jf smOi+ Jj cos Oi (49) 
and the pair interaction 

Uij oc y^sinfy-sinfy) +7j(cos0 / -cos0 J ) +7 l7 cos(0 f - Oj) (50) 

with Jj, Jij, etc. couplings. It is clear that this Hamiltonian (1481) is another example in our 
class, whose partition function can be mapped to the partition function of the AD C/(l) LGT. 

4.4. U(l) LGT with a Monopole Term 

Consider the U(l) LGT with a monopole term. This is described by an additional quadratic 
term supplemented to the standard Wilson action. This system has been studied in connection 
with finding a confinement mechanism based on monopole condensation ll23l l24l l25l . The 
action reads 

S m0 no = J £ (l-cos0 afo (n))+A^|M c (n)|, (51) 

a<b,n c,n 
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where 

M c{n) = J^cdab (®ab(n + d) - ®ab(n)) (52) 

with £ c dab the totally antisymmetric Levi-Civita symbol and n labels the face. The physical 
flux ® a b(n) is related to the face variable ® a b{n) by 

®ab(n) = ® ab (n) - iTt^abin) (53) 

with ^V a b{ n ) the number of Dirac strings passing through the face [|26ll . The last term in (I5TI ) 
has a Fourier series expansion so that this model is included in the complete AD U(l) LGT 

4.5. Mean-Field Theory 

Our proof of the main result yields as a by-product the construction of a mean field theory 
for the AD U{\) LGT This corresponds to an interaction pattern where all pairwise 4-body 
interactions between the edge variables are present. Formally, the action of the mean field 
theory reads 

s m .i.= £ Jijiaco&iei+ej+ek+ei). (54) 

Vz'< j<k<l 

Geometrically, all edges are connected with all edges. This action belongs to the class of 
models whose action we can generate exactly. Therefore our completeness result also includes 
the mean-field theory for the U( 1) LGT. 

5. U(l) Lattice Gauge Theory in a Curved Spacetime 

As a further development and also an important application of the completeness result of 
§12, we study a f/(l) LGT defined in a curved background spacetime and show that it is 
included in our completeness result. Physically, this situation corresponds to photons in a 
curved spacetime, as it happens to cosmological photons that have propagated on a curved 
spacetime during the cosmological evolution. In this section we use the following notation. 
Latin letters a, b etc. are used for flat indices and Greek letters jU, V etc. for curved indices. 
We also adopt the convention that any minus sign before an index can be dragged out as an 
overall (—1) factor because it is abelian, e.g. A_ fl = —A a . The coupling constant is denoted 
by and the lattice spacing by A as in previous sections. 

5.1. Definition of the Lattice in a Curved Spacetime 

In order to formulate the U(l) LGT in presence of an arbitrary background metric g^ v , we 
first define a lattice in a curved spacetime manifold. In the formulation of the U(l) LGT 
in §[2] we started with the Euclidean metric 8 a b = Diag(l, 1,1,1) instead of the Minkowski 
metric r} a b = Diag(— 1, 1,1,1), as the two are related by a Wick rotation. The Wick rotation 
method has been generalized to curve spacetime ll27l . Given a background spacetime metric, 
a smooth family of local Wick rotations can be defined on each (co-)tangent vector space 
of the manifold (thus it is defined on the whole vector bundle). This is done by analytic 



U(l) LGT Connects All Classical Models with Continuous Variables 



17 



continuation of either the local time coordinate or the local metric l|28l |29ll . Therefore, 
Lorentzian manifolds can be Wick-rotated to Riemannian manifolds in general. In our 
formulation below for the curved spacetime, we shall take the metric g^ v to be Riemannian 
(positive definite) rather than Lorentzian for convenience without loss of generality. 

Now let us take a general Riemannian manifold with a given metric g^ v - We could have 
defined a lattice using the coordinate curves (with respect to the natural coordinates), i.e. we 
pick up a set of hypersurfaces defined by x^ 1 = const, such that the manifold is filled up by 
volume cells formed by these hypersurfaces. The intersections of hypersurfaces (coordinate 
curves) induces a 'net' (a graph structure) with 'nods' (vertices) linked by segments of 
coordinate curves (edges). This 'net' formed by coordinate curves might give rise to a lattice 
structure at first sight (once we specify a lattice constant). However, it turns out inconvenient 
to work with the coordinate basisjjj] Therefore, we switch to the non-coordinate basis defined 
by introducing the frame fields e^ such that 

S MV = e^ru, g» v = e£e^ 5 (55) 

where the fields e^ are also referred to as tetrad or vierbein fields which bring the metric to a 
flat one locally. The determinant of the metric field and that of the vierbein are related such 

that g = detfg^y] = ± (det[e^] j with the sign determined by the signature of rj^. Note that 
in our case we have chosen the flat metric to be Euclidean rj a j = 8 a b- Therefore we have 



|det[e£]| = ^de% MV ] = Vg. (56) 

Formally, we have switched from the general coordinate basis 

d 



<9p = -jj-p- 1 for the local tangent space, 

{dx^ 1 } for the local dual space (57) 

with metric g^ v to the orthonormal non-coordinate basis 

{e a = e^d^} for the local tangent space, 

{ (b a = e° dx^ } for the local dual space (58) 

with a flat metric 7]^. With the introduction of the orthonormal non-coordinate basis, now 
we can make use of the integral curves of the basis vectors {e a }, a = 0, 1,2,3 in defining a 
lattice. Note that we do not change the coordination number of the flat or curved lattices. A set 
of integral curves is orthogonal everywhere so that the 'net' formed by these curves defines 
a curved rectangle (rather than an arbitrarily shaped quadrilateral) on each small 'patch' (a 
face). (See Fig. 0) 

Now we could define a lattice by imposing a lattice constant on this 'net' formed by 
integral curves of the non-coordinate basis. It is required that the integral curves are equally 
separated along each direction a, b, etc with a common distance A. The intersections ('nods' 
of the 'net') of this set of integral curves are taken as vertices, so that every segment of integral 

j| The natural coordinates may not be orthogonal and the coordinate differences Ax'' and Ax v (/J ^ v) might 
have different dimensions. It is in general invalid to identify the length of an edge in a lattice with a difference 
in coordinates such as Ax^. 
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Figure 6. The manifold is covered with a net of intersecting integral curves of the non- 
coordinate basis vectors. The curves are equally separated so that each small patch (face) 
is a curved rectangle with equal sides. 

curves in between two 'nods' becomes an edge. Each edge of the lattice has the same length 
A, see Fig. [7J In the figure and in the following we shall use the short notation a = e^Ax^ and 




Figure 7. (a) A small patch (a curved quadrilateral) formed by coordinate curves, (b) A curved 
rectangle with equal sides formed by integral curves of the orthonormal non-coordinate basis. 

\a\ = \b\ = A. 

5.2. The Action of the U(l) LGT in Curved Spacetime 

With a lattice defined, now we can define the edge variables Q a in a similar way as in the flat 
case by 

e a = AWA a (59) 

with 9 a = e« On and A a = e„ A^ . Unlike in the flat case, here we have to be careful and specific 
on the position where exactly the edge variables are defined. It is convenient to put the edge 
variables in the middle of each edge (see Fig. [8]), e.g., take the face with four vertices labeled 
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n, n + a, n + b and n + a + b, the edge variables are written as 

6 a (n + -a), 6 h (n + a + ^b), 0- a (n + l -a + b) , 0_ b (n+-b) (60) 

on the corresponding edges, respectively. 

Next is the face variable & ab which can not be defined by simply taking summation over 
the four edge variables belonging to the face because those vectors are defined at different 
tangent spaces (of different points)^ In order to be consistent with the continuum limit, the 
four edge variables e , e E f for each face have to be parallelly transported to a common point 
before summation. We choose the center of each face so that the expression of 0^ can be 
written in a symmetrical way, see Fig. [8] 

We shall denote by Q e the parallelly transported edge variable Q e from the center of edge 
e to the center of face f (e E f). For example, the parallelly transported 6 a (n + a/2) is denoted 
by a (n + a/2)Q Then the face variable is defined by 

0/=£0 e . (61) 
Take a face as the one depicted in Fig. [8j We have the face variable 




Figure 8. (a) The four edge variables each defined at the center of the corresponding edge, (b) 
The edge variables have to be parallelly transported to the center (or any other common point) 
before summing up. 

®ab = 0a(n + l -a) + e b (n + a + -b) + 0_ £ ,(n + l -a + b) + B_ b {n + l -b). (62) 

The transported Q e is connected with the original Q e by the following relation (taking 
Gb(n + b/2) transported to n + a/2 + b/2 for example) 

e b (n+-b) = e b (n+-b) + -aa}^6 c (n+-b), a = e a ^ (63) 

f The sum Y^eef ®e in the continuum limit A — > approaches A 2 &(F a b + Y^i,A c ) instead of the field tensor F a t 
alone. The extra term involves the structure factor appearing in the commutation relation of the basis vectors 
[i a ,ib] = Y^h^d which never vanishes identically in a non-coordinate basis unless within a finite globally flat 
region. 

+ Edge variables are always transported to the center of the face. Therefore for a given face, the location of the 
transported face is omitted in notation, e.g. In 9 a (n + a/2), the argument n + a/2 only indicates the location of 
the edge variable before transport. 
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with the spin connection co c ab given by 

<, = e c v e^(^+e^ p ). (64) 

Here T p ^ v is the Christoffel connection which is symmetric with respect to the two lower 
indices, i.e. T^ v = Ty^ (torsion free). We emphasize two points here: 

(i) In expressions involving a parallel transport such as (1631) . there is no summation over 
index a if we write e^Ax^ in replacing a in the last term of the right hand side. It 
represents a particular direction rather than a dummy index. Therefore we prefer to write 
directly 

6 b (n+^b) = e b (n+^b) + A(O c ab e c , |a| = |e£A^|=A (65) 

in following parts of the paper (assuming indices a, b etc. not containing a minus sign 
for simplicity). 

(ii) The functions of ab as well as 9 C appearing in the last term are position dependent. 
They could be defined either at the point n + b/2 (as in the expression) or at the point 
n + a/2 + b/2. These different choices have the same continuum limit and are equivalent 
up to ff(A 2 ). 

Therefore, by plugging (1631) into (|62|) . we obtain 

®ab = a (n + X -d) + e b (n + a + hb) + 0- a {n + l -a + b) + 0_ b (n + hb) 

-Aco c ab e c + Acol a e c . (66) 
Let us look at the continuum limit. The sum of the first four terms in (|66l) gives 

eef 

= Ae^(e, v v )-Ae^ v (e^) 

= Ae£e£(d M v - d^) + A(e£^e£ - 4d^ a )0 v . (67) 
On the other hand, using (1641) . we see that the last two terms in (1661) equal 

-A{(o c ab -(o c ba )e c 

= - A (e£(<? M e£ + ejrjj -efC^ + eJry) e c v c 

= - A(e^e£ -^d^l)0y - A(e£ejl^ -e^I^f*, 

= -A(e^e£-e^e fl v )e v . (68) 

In obtaining (l68l) we have used the symmetry properties of the connection with respect to 
the two lower indices. As a consequence of (1671) and (1681) . the face variable in the continuum 
limit approaches 

® ab = J>-A«,-oe c 

eef 

^Ae^(a M v -a v M ) 

= A 2 ^e v b F^ v = A 2 ^F ab . (69) 
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Therefore, with the face variable defined in ( 16"T| ) yielding a proper continuum limit (1691) , 
we defined the action of the U{\) LGT with a background metric as follows: 

^ feF 

This expression resembles the action in © for the flat metric case (with different definitions 
of the faces and face variables). We shall prove that (1701 ) leads to the correct continuum limit. 
As A — > 0, we have each y being small, so that 

s "il j?t= il^F* = l£AVi. (71) 

^ /EF A ^ ab ^ H ab 

Furthermore, we have 

A 4 = |(e«^) A (etAx v ) A (e^P) A (e^Ax^)| 

^ \6D a A 0) b A cD c A cD d \ = |det[e£]|d 4 ;t = ^d 4 * (72) 



being the invariant volume. In obtaining (1721) we have used (1561) with A the wedge (exterior) 
product and 6) a etc. the dual basis. On the other hand, 

F& = F ab F cdn ac n bd = F ab F cd ^ c pg Wz b v 4g Va 

= F^F pa g^g vo = F^ v . (73) 

As a result, by combining ( 1721) and (1731) , we prove that 

A^O 1 



(74) 



Sn°I/dV^ v F<- 

That is, the discrete U(l) LGT action defined in j70l) approaches the pure gauge action 
with a background metric in the continuum limitG 

The expression of the action defined in (1701) with the face variable given in (16TT) involves 
explicitly the lattice constant A because of (1651) and (1661) . It is possible to rewrite the explicit 
expression of the face variable & ab (l66l) in a more similar way to that in the flat case ©. 
Let us look at the last two terms in (l66l) . As shown in (1681) . the Christoffel connections (see 
(l64l) ) cancel out; each remaining term depends on the lattice constant and a derivative of the 
vierbein field in such a way that it is identified as a first order differential of the vierbein field. 
Therefore it can be rewritten as a difference in the vierbein up to the order <^(A 2 ), e.g. 

Ae£(^)e v ~[e> + a)- e/ >)]e v . (75) 

Therefore, the last two terms in (l66l) can be split into four terms (to be more symmetric) each 
involving a difference of two vierbeins, one defined at the center of the face, the other at one 
edge. i.e. 

- A((» c ab e c - (o c ba )e c = -A(^d^ v b -^ b d^)e v 



* 



We should have instead of the field tensor F^ = d^Ay — d v A^, the covariant field tensor J^ v = D^A V —DyA^ 
with the covariant derivative D^A V = d^A v — T^ v Ap appearing in (EOt . However, since the connection is 
symmetric with respect to the two lower indices while the field tensor is antisymmetric, we have F^ v = i£^ v 
identically. 
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~ [e£ (n + l -a + lb) - e£ (n + -a)} 6^ (n + la) 

+ [e^(n+ifl+ife) -$(n + a+lb)}0^n + a+ -b) 
-[4(n + -a + -b)-^(n + -a + b)]e^(n+-a + b) 
- [e£ (n + X -a + hb) - 4 (n + lb)} 6^ (n + lb) . (76) 



Now with this new form (1761) . the face variable can be rewritten in the following form: 

®f = ®ab= J>-A(Q&-a&)e c 

eef 

~ e^(n+la)^(n+la+lb) 
+ 8 ll {n + a+-b)e%{n+-a+-b) 
-8 ll {n+-a + b)eg{n+-a+-b) 
-Op(n + lb)e%(n+la + lb) 

= 52 $edge ' e face- (77) 

eef 

In the last line of (1771) we have used a concise notation to indicate that the face variable 0/ is a 
sum over the edge variables 9 e with each edge variable coupled to the vierbein e« at the center 
of the face. The form (1771) resembles © in the flat case. We also realize that in this form the 
vierbein fields need only to be defined for each face. Moreover, by rearranging terms, this 
form can also be written in the following way 

®ab = [fy(n + la) -e^(n + la + fc)]e£ (n + la + lb) 

+ {e^n + a + lb) - (n + Ib)]^ (n + la+lb). (78) 

The expression (1781) is transparent in taking the continuum limit (reducing to F a b) as well as 
the 'flat' limit (e£ = 5 a M ). 

We emphasize here that a consistent theory of the U(l) LGT with a background metric 
is not unique. (There are two requirements only: 1) The recovery of a U(l) pure gauge 
theory with a metric in the continuum limit; 2) The recovery of a flat U(l) LGT in the limit 
gjiv = l^iv, an d these do not defined the theory uniquely.) We are making a particular choice 
so that the resulting theory takes a form similar to the flat theory and the form of the face 
variable looks more symmetric with respect to the edge variables involved^ 

(j Had we made another choice, e.g. had we preferred the edge variable V in d75l l and 6^ in (IToT l be defined at 
the center of the face rather than on the edge, the face variable could have taken a different form such as 

®ab = % (n + la) -e^{n+ l -a+ l -b)]tf{n + I a ) 

+ [0„ (n + a+ l -b) -d ll {n + \a + I&)]e£ (« + a + l -b) 
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5.3. The Partition Function of the 'Curved' U(l) LGT and the Completeness Result 



We have defined the action of the U (1) LGT coupled to a background gravitational field (a 
background metric). The action is given by (1701) which takes a similar form as in the flat 
case ©. The face variables are defined in (loTI) (by summing over edge variables parallelly 
transported to the center of each face) or in (1771) equivalently (by summing over edge variables 
coupled with the vierbein field at the center of each face). As a result, the partition function 
of the model reads 



■"^curved 




exp 



1 



2^ 2 



1 — cos 




1 dredge 

,e<EE 



exp 



2^ 2 



1 — cos 



y, ^edge ' Sface 



(79) 



The action (TTOl) of the U{\) LGT with a curved background metric satisfies the three 
conditions ©, (flOl) and (fTD) in § 13.11 . Therefore, we conclude that our completeness result 
contains this model, i.e., that the U{\) LGT with a curved background metric can be mapped 
to a4D U(l) LGT with a flat metric. 



6. Generalizations 

The completeness result obtained in § 13.41 can be generalized to include an even broader class 
of classical models. The generalization is achieved by relaxing the third condition (fTTI) in 
§ 13.11 . In other words, we shall consider classical physical systems with dynamical variables 
subject to constraints in this section. A very important set of physical models are classical 
Heisenberg models. For these models the spin variables s are classical but subject to the 
constraint s\ + s\ + ...+ s\ = I, i.e. they live in an internal space which is an sphere. 



6.1. Models with Constraints 

Let us take a classical system satisfying only conditions © and (flOl) . Assume that there are a 



set of M (M < N) independent unsolvable constraints^!! over the variable x/s, such that 

&l{xj} = 0, Z=1,2,...,M. (80) 

The functions &\ are also assumed to have Fourier series expansions. After normalization of 
the variables Xj — > 6j, we write the Hamiltonian as J^{0j} and the constraints as J?{6j}. So 

- [flp (n + )-a + b) - dn(n + \-a + h>)]e£ (» + X -a + b) 

- [0 M (» + ±b) - %(« + ~a + i*)]e£ (n + ±b). 

Also, had we parallelly transported the four edge variables to a common point other than the center of the face, 
the face variable & a /, might have taken a less symmetric form. 

f t If a constraint is solvable, then we could release the constraint by a change of variables such that there will be 
no constraint over this new set of independent variables. 
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that the partition function reads (see (1331) ) 



■^classical eJ Y 



jT flldOH [U 5 J e~^>. (81) 

Now for each constraint j^} we associate an additional variable jq and re-write each delta 
function as a Fourier transform. Then (18TT) becomes 



/ring 



d?Q 

[ 

i 



■ exp | -pj4?{0j} + i£ JQ^i{0,-} | . (82) 

If we make a (large ultraviolet) cut-off on each of the JQ variables such that jq G [—A/, +A/], 
then the integrals f^™ become Now we can normalize these variables by 0/ = jf-JQ. 

Then the partition function (f82l) takes the following form approximately 

^classical ~ o¥ o¥' f \\\^Q^\ e^'^ (83) 



with ^yK' = TJ/ 2^1 an d an effective complex Hamiltonian ,W{Qj, 0/} defined by 

jr{0/,fc} = -js^e^+i^^^e,} (84) 

This Hamiltonian in (1841) with additional variables 0/ has a Fourier series expansion 

<^{0/,0/}= E H {m,,m,}COS I J^mjOj + nnQi 

{ntj,mi} \j,l 

} sin {^ m j e j+ m i < h\ ( 85 ) 

with complex Fourier coefficients H{ OT . and H{ m . jOTi j. Therefore, if we allow for complex 
coupling constants J{ mjm ) = H{ mj: mi} and J{ mjjn ,} = H {m , m} , the partition function Classical 
in (l8~3l) is of a t/(l) LGT type. As a result, if we assume that the cut-off (a regularization) 
in the delta function representation above works universally, then the AD C7( 1) LGT with 
general complex coupling constants is a complete model for classical statistical models with 
continuous variables subject to constraints. 



6.2. Models with discrete variables 

Another possible generalization is to consider classical models with discrete degrees of 
freedom, such as the Ising model, the Potts model, the vertex models, and the 7L q LGT where 
the dynamical variables take values only in a discrete set. These models can be considered a 
special case of continuous models with constraints. 

Let us again consider a classical system satisfying only conditions © and (flOl) . Each of 
the variables xj can assume only finite discrete values 

Xj=X h , aj<X h <b h /; = 1,2, (86) 
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The possible values X[. and the number of possible values Lj may depend on j. It is equivalent 
to express these constraints by the following expression with 5-functions 

Li 



n 



.0=1 



(87) 



such that for an arbitrary function of a set of x/s, we have 



n n 



/{*;}= E (88) 

In the last summation of < [8~8l) we are summing over all possible configurations (all possible 
combination of Xj values). Moreover, after normalization of the variables Xj — > Gj, J^{xj} — > 
Jf{6j}, the delta functions become 



2% 



b j~ a j 



8 6 



2% 



aj + bj 



We define 



0, 



2% 



b i~ a i 



x, 



aj + bj 



(89) 



(90) 



for notational convenience. Consequently, the partition function of a classical model with 
discrete degrees of freedom can be expressed as 

•^discrete 6 ^ ' (91) 



all {lj} 



fk n 



Now, with these particular forms of 5-functions appearing in (|9TT) , it is not necessary to go to 
the region of a complex Hamiltonian. Let us represent each 5-function by a limit 

8 (ej-® h ) = \hn8 Ej (0,-0/.) (92) 



with each 8 £j (0j — 0/.) being a positive-definite function on [ — tt, tt] depending on the 
parameter £,-. (This function could be a Gaussian or a 'square' function.) If we make a 
cut-off approximation on these parameters, then indiscrete can be rewritten as 

•^discrete (93) 



(2w) 



1^,(0,-0/,) 



The expression on the exponential of ( |93l can be interpreted as an effective (real) Hamiltonian 
which allows a Fourier series expansion over 0,'s with real coefficients. Therefore we 
conclude that U(l) LGT is complete also for certain classical discrete models, again assuming 
that the cut-off (regularization) (|92l ) works universally. 
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7. Conclusion 

We have proven that any classical partition function depending on continuous variables subject 
to conditions ©, (flOl) and (fTTI) can be approximated (to an arbitrary precision) by the partition 
function of a AD U{\) LGT. In the proof we first introduced a quantum representation of the 
U{\) LGT partition function. Then through merging and deletion of gauge field variables 
and proper choices of local coupling constants, a mapping from a AD U{\) LGT partition 
function to a more general partition function is established. In this sense the AD U(l) 
LGT is a complete model for a large class of classical models. The completeness result is 
also generalized to include continuous models with constraints (if we are allowing complex 
coupling constants) and discrete models. As a further development and important application 
of the completeness result, we have developed a consistent theory of the U{\) LGT coupled to 
a background metric. The action is defined in a form (|7Q|) very close to that of the model in flat 
spacetime. Our completeness result holds for this model such that its partition function can 
be mapped (approximated to an arbitrary precision) to the partition function of a U{\) LGT 
in flat spacetime. This is the first time that a completeness result is established for continuous 
statistical models. 

We believe that our completeness result cannot be proven with a 3D U(l) LGT. Another 
open question is whether a similar result can be found for non-Abelian LGTs (e.g. for SU (2) 
or SU (3) LGTs). We envisage that these theories may require a new approach since a direct 
generalization of our construction seems not to be possible. 
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Appendix A. Construction of Many-Body Interactions 

To generate the functions of the set (|T4|) we will proceed similarly as in [J3l SJ, namely we 
will make use of the merge and the deletion rule, and of the gauge fixing of edges. The latter 
is a procedure for fixing the values of the variables in the lattice, which results in a theory 
that is physically equivalent to the original one. This can be carried out by virtue of the 
gauge symmetry of the model. The only restriction in this procedure is that the edges whose 
variables are fixed by the gauge cannot form a closed loop j[22|. 

To construct the many-body interactions, we first need to "propagate" the variables 
inside the lattice in order to bring them close together to interact. This propagation is 
achieved with the following construction (see Fig. IA1I) . On a cube (on the left of Fig. IA1I) . 
we merge the face at the front, at the bottom, and at the back to generate the interaction 
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Jfcos(oci — OCi + 0i + 62). Then, we let this coupling strength go to infinity, Jf — > °°, which 
imposes the constraint 61 + 62 = 0. The same process is repeated for the cube on the right of 
Fig. IA11 where the constraint reads 62 + 63 = 0. Thus, we have set 63 = 6\, that is, 6\ has 
"propagated" two sites to the right. Note that if 6\ were propagated an odd number of times, 
the resulting variable would equal —6\. This can be circumvented by letting 6\ participate 
in the final interaction with the opposite sign (as explained below). In order to see how to 
turn the propagation path, we refer the reader to the explanations on Fig. 11 of JH, since the 
construction is analogous. 




Figure Al. Propagation of the variable Q\ across the 4D square lattice (the figure shows only 
a 3D projection of this space). In all figures, red edges denote edges whose variable has been 
fixed by the gauge. By means of the merge rule, Q\ is propagated into 63, since O2 — — Q\ and 

63 = —02- 

Now we focus on the replication of the (classical, continuous) variables, that is, on the 
generation of several copies of a given variable. This is achieved by applying the propagation 
procedure explained above into the fourth dimension, as shown in Fig. IA2I The interaction in 
the yellow cube is of the form Jf cos (62 + OC3 + 64 — OC3), on which we let Jf — > °° and thereby 
impose 62 = —64 (i.e. we apply the merge rule on this cube as well). The rest works exactly as 
the propagation explained in Fig. lAll We note that the reason for using a fourth dimension in 
the replication of edge variables is that all schemes we have found in three dimensions involve 
closed loops of variables fixed by the gauge 0. 

Next we show how to generate interactions of the type (fl4l) . The generation of the 
interactions cos( Y*f=\ &i) f° r the specific case K = 5 is illustrated in Fig. lA3l We first propagate 
each of the variables 6\,. . .,65 inside the lattice in order to distribute them on the edges of 
the rectangular prism shown in Fig. IA31 Then we merge all faces on the exterior surface of 
this prism into one large, blue face. As shown above, this face only depends on the spins at 
its boundaries, that is, it has the form 

//COS(01 - «2 + 02 - OC3 + OL4 - «4 + «3 + «2 - «7 

+ a 6 + a 5 + 05-a5 + 04-a 6 + 03 + a 7 ). (A.l) 

Since the dependence on every auxiliary variable a ; cancels, it takes the desired form 
^/ C0S (XjLi 0;)- The generalization to any K is straightforward. For odd K, we construct a 
longer or shorter prism than that of Fig. IA31 arranging (K — l)/2 variables in front of the 
remaining (K — l)/2 (as 6\ and 62 are arranged in front of 63 and 64). The remaining variable 
would be unpaired, and the vertical edge on the corner of the prism fixed by the gauge (just 
as 65). For even K the construction is simpler. For K — 4 we arrange the spins 0\, . . . , 64 as 
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Figure A2. Yellow faces are in the direction of the fourth dimension and have the same 
meaning as blue faces (i.e. merged faces). Replication of the variable Q\ into 03, 65 and 06. 
This replication is essentially a propagation (as the one of Fig. IAU in the fourth dimension, 
i.e. the variable lives now in another 3D space. 



in Fig. IA3[ and we merge the large, blue face over the face joining the red u-shapes of 62 
and 04 on the right. For any other even K, we similarly arrange K/2 variables in front of 
the remaining ones. Finally, note that the generation of interactions of the type cos(20i) is 
achieved by first replicating the variable 9\, and then letting its two copies participate in a 
two-body interaction as explained above. 
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Figure A3. Bold, black edges contain variables that participate in the final interaction 
01 , . . . , 05. Red edges (except for the one marked with a) stand for edges whose variables 
have been fixed to zero using the gauge symmetry, and blue faces stand for merged faces, as 
in Fig. [5] If a = 0, as the other spins, the large blue face contains the five-body interaction 
7f cos(0! + 2 + 03 + 04 + 05), whereas if a = -n/2, it depends on 7/sin(0i + 2 + 03 + 04 + 

05). 



Now we only need to show how to generate interactions involving variables with different 
signs. The generation of 7ycos(0i + 02 + 03 — 04 — 0s) is shown in Fig. IA4I The variables 
which have the same relative sign are arranged as explained for the cosine of the sum of 
variables (see Fig. IA3I) . The new element here is that the two sets of variables which have 
the opposite relative sign must arranged perpendicularly to each other. Then a large, blue 
face is merged over the external faces, where the desired interaction takes place. One can 
also verify that the dependency on the auxiliary variables 0C ; cancels out, and that variables 
perpendicular to each other have the opposite sign. The generalization to an interaction 
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JfCOs(Y,f=i 0/ — X^Li Qj) 1S a ^ so straightforward. One only has to arrange the first set of 
variables as explained for the case Jfcosfcfli 9i), with K\ either odd or even. The other 
set of variables is arranged also as explained above (with K2 being odd or even), and set 
perpendicular to the first set. 




Figure A4. Five-body interaction Jf cos ( Q\ + 62 — 63 — 64 — 9s ) . The meaning of the symbols 
is the same as in Fig. IA3I 



Finally, we point out that the sine functions are generated by making use of the relation 
sin(y) = cos (7— n/2). This phase amounts to gauge fixing one of the spins to —n/2 instead of 
0, as the rest. For example, to generate sin($i + 62 + 63 + 04 + 65) we construct the interaction 
of Fig. |A3l and we fix o = —n/2. 

This shows how one can generate all interactions of the set (fl4l) starting from a AD U(l) 
LGT. The construction also shows that each Fourier basis function (i.e. each term in the set 
([14)) ) requires a polynomial enlargement in the number of variables of the AD U(l) LGT. 
We shall return to this fact in § 13.41 . where we specify the overhead in the system size of the 
complete model as a function of the features of the target model. 



Appendix B. Accuracy of the Finite Fourier Series 

In the completeness result of § 13.41 . we have made a truncation in the Fourier series basis (see 
(|3TI) and (|32|) ) so as to approximate a general Hamiltonian as an expansion. Here comes a 
question of accuracy. That is, given a truncation — M < mj < M, Vj in the Fourier modes, 
how close is the following finite Fourier series 



H 



L<!} ex P i!>./% CB.l) 

{mj} I J J 

to the original Hamiltonian function H^ K \{9j}) of a A"-body interaction term? According 
to ll30l . for a smooth enough (usually at least differentiable) single variable function f(6) 6 
C a [-n, ft] with a > 0, we have 

\^M[f]-f(e)\<A(f) < B - 2) 

where [f] is a finite Fourier series of f(6) with 2M+ 1 terms (truncated at the M th mode as 
in our case), A(f) is a finite factor depending on the function form of f(6) on the domain only, 
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and the finite sum is over some polynomials P(oc) of a. Therefore we see that in order to have 
an accuracy ~ for the finite Fourier series of a single variable, we need a polynomially 
large truncation Poly{^Y) in the Fourier modes. We could say that we have a polynomial 
accuracy in this case. Now with a K-body interaction term which is a function of K variables, 
in order to have an accuracy 

-tf w ({e,})|< Jp, (B.3) 

we shall need J% ~ [Poly(jV)\ terms in the finite Fourier series, i.e. the K th power of some 
polynomial of jV . Therefore it is still efficient if the model has only few-body interactions 
{e.g. nearest neighbor interactions only) or if K scales polynomially with the system size. 




References 



Pathria R K 1996 Statistical Mechanics (Elsevier, Second edition) 

Creutz M 1983 Quarks, gluons and lattices (Cambridge University Press, Cambridge, England) 
De las Cuevas G, Diir W, Briegel H J and Martin-Delgado M A 2009 Phys. Rev. Lett. 102 230502 {Preprint 
larXiv:08 12.35831 1 

De las Cuevas G, Diir W, Briegel H J and Martin-Delgado M A 2010 New J. Phys. 12 043014 (Preprint 
larXiv:091 00961) 

Van den Nest M, Diir W and Briegel H J 2008 Phys. Rev. Lett. 100 1 10501 (Preprmf larXiv:0708.2275l l 
De las Cuevas G, Diir W, Van den Nest M and Briegel H J 2009 J. Stat. Mech. P07001 (Preprint 
larXiv:08122368l) 

Van den Nest M, Diir W, Raussendorf R and Briegel H J 2009 Phys. Rev. A 80 052334 (Preprint 
larXiv:0805.1214l) 

De las Cuevas G, Diir W, Briegel H J, Van den Nest M and Martin-Delgado M A 2010 Quantum Algorithms 

for Classical Lattice Models (in preparation) 
Hubener R, Van den Nest M, Diir W and Briegel H J 2009 J. Math. Phys. 50 083303 (Preprint 

larXiv:0812.2127l) 

Van den Nest M, Diir W and Briegel H J 2007 Phys. Rev. Lett. 98 117207 (Preprint 



arXiv:quant-ph/06 10157) 



Panero M 2005 JHEP 0505 066 (Prepnwf |arXiv :hep-lat/0503024); 2005 Nucl. Phys. Proc. Suppl. 140 665 

(Pre/?n>if |arXiv:hep-lat/0408002| l 
Wilson K G 1974 Phys. Rev. D 10 2445 
Polyakov A M 1975 Phys. Lett. B 59(1) 79 
Polyakov A M 1975 Phys. Lett. B 59(1) 82 
Wegner F 1971 J. Math. Phys. 12 2259 
Kogut J B 1979 Rev. Mod. Phys. 51 659 
Gross D J and Wilczek F 1973 Phys. Rev. Lett. 30 1343 
"t Hooft G 1981 Nucl. Phys. B 190 455 

Cahill K and Reeder R 1986 Phys. Lett. B 168 381; 1986 Jour. Stat. Phys. 43 1043 
Palumbo F 1990 Phys. Lett. B 244 55 

Patrascioiu A, Seiler E and Stamatescu I O 1981 Phys. Lett. B 107 364; Seiler E, Stamatescu I O and 

Zwanziger D 1984 Nucl. Phys. B 239 177. 
Creutz M 1977 Phys. Rev. D 15 1 128 
Barber J S and Shrock R E 1985 Nucl. Phys. B 257 5 15 

Damm G and Kerler W 1998 Nucl. Phys. Proc. Suppl. 63 703 (Preprmf [arXiv:hep-lat/9709061| l 
DeGrand T A and Toussaint D 1980 Phys. Rev. D 22 2478 
Dirac P A M 1931 Proceedings of the Royal Society A 133 60 



U(l) LGT Connects All Classical Models with Continuous Variables 



31 



[27] Birrel N D and Davies P C W 1984 Quantum Fields in Curved Space (Cambridge University Press, New 
York, NY) 

[28] Liu C-H 1997 (Pre/?ri»f [arXiv:hep-th/9707"T9^ l 

[29] Moretti V 2000 Commun. Math. Phys. 212 165 (Preprircf |arXiv:gr-qc/9908T)6l^ 
[30] Schlag W Notes on Harmonic Analysis http://math.uchicago.edu/^schlag/book.pdf 



